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Abstract. We propose to study the L 2 -norm distance between classical and quantum 
phase space distributions, where for the latter we choose the Wigner function, as a 
global phase space indicator of quantum-classical correspondence. For example, this 
quantity should provide a key to understand the correspondence between quantum 
and classical Loschmidt echoes. We concentrate on fully chaotic systems with compact 
■ (finite) classical phase space. By means of numerical simulations and heuristic 

arguments we find that the quantum-classical fidelity stays at one up to Ehrenfest- 
type time scale, which is proportional to the logarithm of effective Planck constant, 
and decays exponentially with a maximal classical Lyapunov exponent, after that time. 
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1. Introduction 

In this paper we want to address few questions concerning the correspondence of 
quantum and classical evolution on the classical phase space for quantum pure states 
of classically (strongly and weakly) chaotic systems. Our description is restricted to the 
compact or effectively compact phase space. By effectively, we mean that the classically 
available space is compact and quantum description can be effectively done in some 
finite dimensional set of basis functions. 

We will use the Weyl- Wigner phase-space representations of quantum mechanical 
states jl] and in this way compare the classical and the quantum evolution, both 
as evolution of functions defined over the classical phase space. The analysis will 
be performed by considering coherent packets as initial quantum states, such that 
the corresponding Wigner functions coincides with the initial classical phase space 
distributions identical to Gaussian packets. It is a common knowledge in the realm 
of quantum chaos that the phase space correspondence between classical and quantum 
mechanics drops down on the scale of Ehrenfest time £e = — log7i/A, where h is an 
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effective Planck constant h = ^physical/A A is a typical action of the system and A is the 
maximal Lyapunov exponent. For a compact phase space the Ehrenfest time is usually 
defined as £e = log N/X, where N is the Hilbert space dimension. Here we would like to 
investigate in more detail the break down of the correspondence and the corresponding 
break time scale. 

Let us systematically define the problem. First of all, we start with a classical 
Hamiltonian system defined over a given (i.e. sympletic ) phase space X, for example 
X = T*R™ or X = T 2 . The trajectory of this system at the continuous (discrete) time 
t G R + (t G N U {0}), initially in x G X, is given by 0'(x) G X. The corresponding 
quantum system is defined over the Hilbert space TC and is described by a state \ip) G 7i. 
The motion of the quantum state is given by the unitary evolution operator U l , with 
t G R + and t G N U {0}, respectively. 

To represent the quantum operators, A : H — * H, as functions over the phase space 
X we use the Wigner representation A — > H / [A](x) which share the same main properties 
in all classical phase-spaces, even if its specific form strongly depends on the geometry 
of X. We will recall the exact definition of these quantities on Euclidean X = M? d and 
toroidal phase space X = T 2 . We are only interested in the phase space representation 
of density operators of pure states p = called the Wigner functions and denoted 

by W^(x) = W[p\(x) G E and normalized as J x dxW^(x) = 1. By using the Wigner 
function one can represent quantum states and their dynamics in terms of phase-space 
functions. Similarly we can describe the state of a classical system on the phase space 
in terms of probability density functions 

We study the classical- quantum correspondence (CQC) by observing the deviation 
between the Wigner function of the quantum system and the corresponding probability 
density of the classical system evolving in time. The initial state ^(O)) of the quantum 
system is a coherent state of minimal width with the Wigner function W t=0 (x) depending 
on the geometry of X, but is basically similar to a symmetric Gaussian probability 
distribution! or is converging to its form in semi-classical limit h — > with error 
approximately of order exp(— 10(^ _1 )|). The last estimate is refined in the following 
sections for the toroidal geometry using the reference [3]. The state of a classical system 
at t = 0, denoted by p t=0 (x), is a Gaussian probability distribution corresponding to 
W t=0 (x). The phase-space functions of the initial quantum and classical states are either 
identical, or converge to each other in the semi-classical limit H — > 0. Our aim is to 
understand the behavior of the L 2 norm of the difference between the time evolving 
Wigner function of the quantum system W t (x) = Wut\ip(p)){ x ) an d the time evolving 
probability density of the classical system p t (x) = p°(0~'(x)), defined as 



The evolutions of the Wigner function W t (x) and of the classical density p t (x) are 
both unitary so L 2 norm of both functions is conserved: 1 1 /?* ( ) 1 1 2 = II P°{x) \\ \ and 

% On compact phase spaces it is sometimes necessary to periodicize the Gaussian distribution to obtain 
maximal similarity with the Wigner function of a coherent state. 




(1) 
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II ^l 3 -) II 2 = Il^°( x )ll2- The only time evolving quantity in the expression (JTJ) is the 
overlap of the Wigner function W l and the classical density p* named the quantum- 
classical fidelity (QCF): 

F(t) = [ dxp t {x)W t {x). (2) 
Jx 

The QCF (j2J is the measure for the quality of the CQC The F(0) is the maximal value 
of F(t) on Euclidean phase space, or is maximal in the limit h — > for compact phase 
spaces. Therefore, in order to study CQC we observe relative QCF defined as 

G(t) = |||<l + exp(-|0(r 1 )|). (3) 

On the Euclidean phase space we can skip the h corrections on the right hand side of (jHJ • 
The relative QCF G(t) is the main object through which we study CQC in the following 
sections. The decrease of G(t) from the value 1 is understood as notable violation of the 
CQC. 

An important motivation for the study of QCF comes from the studies of classical 
[I] and quantum Loschmidt echoes and their correspondence. One can immediately 
show that a good control of QCF as defined above (J2J) gives in turn a good control of 
the correspondence between the two Loschmidt echoes. In order to define the Loschmitd 
echoes, or fidelities, we have to follow two slightly different time evolutions run by two 
slightly different Hamiltonians. Let the classical Liouville densities and quantum Wigner 
functions of the unperturbed and perturbed evolutions at time t be denoted as Pq(x), 
Wq(x), and pl(x), W*(x), respectively, where e characterizes the size of perturbation. 
Then the classical Loschmit echo (CLE), or classical fidelity, and the quantum Loschmidt 
echo (QLE), or quantum fidelity, are defined respectively as 

F? hE (t)= j dxpKxUix), (4) 

F QLE ^ = Jdx Wl (X) Wt (X). (5) 

For strongly chaotic systems, within the Ehrenfest time scales, QLE has been established 
to exhibit an exponential decay with the rate given by a classical Lyapunov exponent jSj. 
On the other hand, again for sufficiently short time scales, CLE also exhibits Lyapunov 
decay as has been derived in p]. These results may suggest a simple conclusion that 
the Lyapunov regime of QLE is a purely classical phenomenon. Yet, there is a subtle 
issue of the precise time-scale (perhaps a fraction of the Ehrenfest time) up to which the 
connection may be established. And this connection can be given through the concept 
of QCF as shown below for the Euclidean geometry X = M. 2d with d degrees of freedom. 
In order to do that we shall make use of an elementary inequality between arbitrary 
quadruple of Hilbert space vectors u,v,r,s 

\\u — v \\2 — 1 1 J" — s || 2 | < || tt — r|| 2 + ||l> — S || 2 (6) 

The inequality can be proven straightforwardly by means of standard triangle 
inequalities. Now, putting u := Pq, v := p', r := Wq, s := W*, and assuming that 
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Wigner-functions and classical densities are all square normalized to a unit Planck cell 
H M I|2 = IMI2 = Iklb = IHI2 = (2nh)^ d due to purity of the state, we obtain an upper 
bound estimate on the difference between quantum and classical fidelity in terms the 
sum of QCF of unperturbed and perturbed dynamics 

1 - F e CLE (t) - ^l-F? LE (t) < ^Jl-F® CF (t) + ^l-F e QCF (t). (7) 

For example, if within a certain time t, QCF remains close to 1 for both evolutions, then 
within the same time QLE closely follows CLE. 

In this paper we discuss CQC on Euclidean and toroidal geometry, the latter being 
an example of a compact phase space. Toroidal geometry is particularly useful, because 
in this case we know plenty of simple chaotic models with simple exact quantizations and 
additionally it supports a smooth and a discrete Wigner function (WF) formulations. We 
study CQC numerically through the QCF using the Sawtooth map and the Perturbed 
cat map as canonical examples for the linear-discontinuous, and nonlinear-smooth maps, 
respectively. The smooth WF on the torus can be interpreted as a model for WF defined 
on Euclidean phase space on which me make some heuristic predictions on the decay of 
the QCF. 

Some related comparative studies of classical and quantum phase space 
distributions in the relation to Lyapunov chaos can be found in Refs. |H1 E] • 

The paper is organized as follows. In section 2 we outline a heuristic derivation of 
the Lyapunov decay of QCF in the echo picture. In section 3 we give more rigorous 
definitions of QCF for the case of compact two-dimensional toroidal phase space and 
two possible distinct definitions of the Wigner function, and later on present numerical 
calculations for two specific models, namely Sawtooth map and Perturbed cat map. In 
section 4 we summarize our results and formulate conjectures for future rigorous work 



2. Euclidean phase space and heuristics 

The Euclidean (symplectic) classical phase space is perhaps physically the most relevant 
geometry and is here denoted by X = Q x P = M. 2d , where Q = {q G M. d } is the 
configuration and P = {p G M. d } is the momentum space. By using some heuristic 
derivation in this geometry we qualitatively explain the numerically observed behaviour 
of QCF (J2J) for the situation, when classical and quantum system are initially described 
by a Gaussian packet for a generic chaotic system with effectively finite available phase 
space. The finiteness of the phase space assures, due to Poincare recurrence theorem fTU] . 
that quantum wave-packet stretched due to chaotic dynamics eventually interferes with 
itself causing a clear violation of the CQC. The point in time when this happens we call 
initial break of CQC and we are mainly interested in the development of QCF after that 
point in time. Without loss of generality we can restrict ourselves to time independent 
Hamiltonian systems. We further restrict our discussion to a simply connected chaotic 
component of phase space (e.g. the energy surface). 
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The classical and quantum system are defined in terms of the Hamilton function 
H(x), x £ X and the Hamilton operator H : TC — ► TC over Hilbert space TC, respectively. 
We shall be working in Weyl-Wigner quantization. The classical dynamics can be 
described solely in terms of characteristics of the flow, namely the classical trajectories. 
We write a symplectic map 0* : X — > X which represents a trajectory starting in 
the point x G X as function The quantum dynamics does not support such a 

treatment and we usually refer to this property as non-locality of quantum mechanics. 
The classical propagation of the phase-space probability density p(x) is defined by the 
Liouville equation 

p = Cp, e tc p(x)=p(r\x))=p t (x), £ = {H,.}, (8) 

Here we use the Liouville operator C, which we express using the bi-directional 
differential J as 

c = h7, 7=$>^f, j= K * 

ij L 

where superscripts L (left) and R (right) indicate the direction in which the differentials 
act and matrix J is the symplectic unit. The Wigner-Weyl formalism in this geometry 
is well known |L. The Wigner function W^(x) of the quantum state \tp) is defined by 

w *( q >ri = jdny J dv( ? pv/n &\<i + v / 2 )(<i- v / 2 \^- ( 10 ) 

Similarly to classical propagation of the probability density, we can define the dynamics 
of the Wigner function W t (x) using the Weyl-Liouville operator C w introduced by Moyal 
1J as 

W l = C w W\ WHx) = e tc ™W(x), C w = H — sin(aj) , (11) 

a 

where we introduce a constant a = h/2. The operators C and £ w are generators of 
unitary dynamics on L 2 (X) and they do not commute in general [C,C W ] ^ 0, which is 
essential for the properties of quantum mechanics. The coherent state with a deformation 
cr in position direction and centered at the phase space point ((fcbPo) G X is written as 



(?l(lb,Pb>= (j-) e-arto-»r+-^ . (12) 

The product of uncertainties in position and momenta reads AqAp = h/2. The Wigner 
function corresponding to the coherent state (fT2"j) is 

WW*) = (^p e " f (9 " 90)2 "^ (P " P0)2 . ( 13 ) 

which is completely positive with a unit phase space integral. The functions of the form 
(JT3J) are used as initial conditions in the evolution of classical and quantum system, 
represented by the Wigner functions W t (x) and the classical probability densities p t (x): 

W t = (x)=p t =°(x) = W {qo , Po) (x). (14) 
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By considering the time evolution of the Wigner functions and the corresponding 
classical density function the QCF (J2J) takes in this geometry the following explicit 
form 

F(t)= [ AxW Q {x)[e- tc e tc ™W°]{x) = [ dx z°{x)z t {x) , (15) 
Jx Jx 

which we simplify by introducing a convenient observable z t {x) = e~ tc W t named 

classically echoed Wigner function. The dynamics of z t (x) is defined by equations 

z = H(t)z, z t = e-'V £w z°, TZ{t) = e- tc (C w - C)e tc , (16) 

with the initial condition z t=0 = W t=0 . The generator of evolution lZ(t) is just the 
difference between the quantum and classical generators, C w and C, in the classical 
interaction picture and is explicitly written as 

*W = E E Mdt +1 Hmx)) e^df^ . (17) 

n=l {Zn+ l) - ij=l 

where we assume smallness of the semi-classical parameter a <C 1. In the strong chaotic 
case with sufficiently simple topology of the flow the leading contributions to TZ(t) can be 
estimated by considering that pre-factors (<9 ? m if )(0*(o;)) behave as stochastic variables in 
time and that hyperbolicity of the flow causes the derivatives in the interaction picture 
to grow exponentially with the asymptotic rate given by maximal Lyapunov exponent 
Amax as lle^c^e^H ~ exp(A max £) [Zj. Note that the observation time should be long 
enough so that the unstable direction clearly develops, meaning A max t ^> 1. On the other 
hand, it is meaningless to treat H(t) (fTTj) for times much larger then i E = | log/z|/A max , 
because for t ^> t E the Wigner function already spreads over the whole space. Taking 
these facts into account we find that QCF F(t) can be qualitatively estimated by 

||1 _ e |loga|-A maxt ^ l <<; Kgxt < | logQ[ | (lg) 

We conclude from the formula ()18j) that the QCF decays exponentially with the rate 
given by the maximal Lyapunov exponent. In classically ergodic systems the probability 
density converges to the invariant density denoted by /i(x): 

lim p t (x) = fx(x) , e tC fi(x) = fi(x) , 

t— >oo 

and therefore one expects that QCF, for times t ^> t^, converges to the classical ergodic 
average (•) of the Wigner function 



lim F(t) = f dxW\x)fx(x) = (W t ) ( 



The ergodic average of the Wigner function can be expected to be roughly constant 
in time. This means the QCF drops to a plateau. However, for systems with relatively 
small Lyapunov exponents, namely if i E is large and comparable with the Heisenberg 
time tu ~ 1/h, QCF could experience another regime of relaxation towards the plateau 
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related to decay of classical temporal correlation functions. This is easily understood, 
if we take into account two things. First, the classical density p t (x) converges towards 
the invariant measure with the rate determined by the leading Ruelle resonance v [T3] § 
with the corresponding eigenfunction (f) v { x\ clS 

p\x) « p{x) + <j> v {x) ((/> v , p°) e~ ut , 

and second, that at time t ~ in the wave-function of a classically chaotic system is 
almost a random function with fluctuations on the scale of h~ and similarly for the 
corresponding Wigner function ^3]. Putting this two facts together, we can estimate 
the asymptotic decay of QCF down to the plateau by 

F(t) - (W) e « e~ vt (<p u , p°> (W*, <p u ) . (19) 

By assuming that the Wigner function lives in effectively finite Hilbert space 7i and by 
introducing the projector P : L 2 (X) —>■ 7i from full Hilbert space to the effective finite 
one, we can approximate the matrix element 

l<wj,^>| < ||p<M 2 . 

3. QCF on the torus T 2 

3.1. Quantum mechanics over the torus T 2 

The mathematical description of quantum mechanics over a compact phase space has 
a long history. Here we restrict ourselves to recall the main definitions and properties 
that are need in the exploration of quantum-classical fidelity, refereing to ^3] and [21] 
for details and references therein. 

The most important consequences of doing quantum mechanics on a compact phase 
space such as T 2 = [0, l] 2 are that the inverse of the Planck's constant h~ can only take 
positive integer values and moreover the corresponding Hilbert space of quantum states 
turns out to be finite dimensional with dimension exactly N = h~ x . This is an immediate 
consequence of periodicity of states in both configuration and momentum representation 
and it will have, as we will see, important consequence on the structure and properties of 
both coherent states and Wigner functions. Just for the purpose of fixing the notations, 
let us recall briefly the construction of the Hilbert space of states and the implementation 
of the dynamics. 

We start by defining the quantum Hilbert space associated to the torus phase space 
T 2 . Physically we must select wave functions satisfying the same Bloch conditions: 

i>(q + L)= e 2M ^{q), ^{p + L) = e~ 2M i,{p) , 

§ In some dynamical systems the Ruelle resonances are well established term i.e. Axiom A, geodesic 
flow on a manifold of negative curvature, whereas in many others, even with exponential decay of 
correlations, the situation is not yet clear. 
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where we have used the ft-Fourier transform 

i>(p) = (F*0)(p) := [°° ^ q )e-^ h -^= . (20) 

For this short presentation of the quantum kinematics we set L — 1 and use periodic 
boundary conditions Q\ = 62 = to keep the notations simple. Later on, in order to 
simplify the quantum propagator associated to the sawtooth map, we will consider L 
even and it is important to stress that all the considerations we make here can be 
naturally extended to the more general case of a torus with arbitrary integer sides. 

Note that we can characterize periodic wave functions using the Heisenberg 
translations: for any h G (0, 1], we consider the quantum translations (elements of the 
Heisenberg group) T v = Q 1 ( v ^- v ip)/ h ^ v G M 2 , acting on L 2 (IR), obtained by the canonical 
quantization of the linear translations in phase space, with infinitesimal generators q 
and p. By extension they also act on a larger set of wave functions, not necessarily 
in L 2 (R), such as periodic functions or infinite linear combinations of delta functions 
over a discrete periodic lattice. In this distributions space «S'(M), we define the space of 
invariant wave functions 

Tin = G <S'(R) : T (li0) ^ = T (0il) ^ = ^} • 

These are distributions ip(q) which are Z-periodic, and such that their ft- Fourier 
transform (jSOj) is also Z-periodic. One easily shows that this space is nontrivial iff 
{2nh)~ l = N G N, which is assumed from now on, and is denoted by TCn = 'Hn- It forms 
a ^-dimensional vector space of distributions admitting a "position representation" 



N-1 



N-1 r . 1 N-1 

1- 



3=0 



*i \Qi) > (21) 



3=0 

where each coefficient tpj G C. Here we have denoted by {lOj)}^) 1 the canonical 
("position") basis for Hn- In the same way as we produce generic quantum states 
over the torus by "periodizing" quantum states over the line, we will construct the basic 
objects of our study, such as quantum propagators, coherent states, Wigner functions, 
etc... by again "projecting" down to the torus the similar object on the line or on the 
plane. In particular let us define the "projector" on L 2 (IR) 

A* = E (-i) 7 "™ f m =(J2 ( E 4*,o)) . (22) 

meZ 2 m 2 eZ mieZ 

Then TL^ = Pj^ («S'(R)) and we will use in particular this projector for defining coherent 
states on T 2 . Given now a classical smooth observable / G C°°(T 2 ), we want to describe 
how to associate to / a suitable quantum operator / on TC^. On the "usual" planar R 2 
case, there are two main quantization methods that can be both adapted to the torus: the 
Weyl and the (positivity preserving) anti-Wick quantization based on coherent states. 
For the purposes of this paper, we will here briefly describe the first one, refereing to 
[20] for what concern the second one. Note first that any smooth function f(q,p) on the 
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(unit) torus can be considered as a double periodic function over the plane R 2 and it 
can be expressed via the discrete Fourier transform: 

f(q,p)=Y.f( k y Mqk2 ~ Pkl) - (23) 
kez 2 

The Weyl quantization of this function is the following operator: 

f=J2 />) T ( k ), where T ( k ) = T k/N . (24) 
kez 2 

It is immediate to see that fH-N ^ 'Hn- The quantization map /—>•/, due to the finite 
dimension of the quantum Hilbert space, turns out to be a non invertible transformation. 
In fact, because of Heisenberg relations T(k + mN) = (_i ) fe 2mi-m 2 fci T(k), we have the 
periodic relation T(k + 2miV) = T(k). Namely, letting Z 2iV = {0, 1, . . . , 2N - 1} and 
assuming / smooth: 

/= E (E/>+ 2£iV )) T ( k )- ( 25 ) 



This, together with (J23)) and fast decay of the Fourier coefficients, immediately implies 
f = g for any given two smooth classical observable f,g G C°°(T 2 ) that agree on the 
2iV-grid classical lattice in phase space : f{^i^h) = iK-iv"' w)> Vm, n G Z 2 jv- As it 
is discussed for example in [22], this fact will also effect the properties of the Wigner 
function associated to a given quantum state|| \I/ G 7^at, that nevertheless will share 
the fundamental phase-space distribution property: 

W>> = N [ W*(q,p)f(q,p)dqdp. (26) 

Before coming back to some more details concerning Wigner functions and coherent 
states on the torus, let us now briefly recall the main definition and properties of the 
quantum dynamics over T 2 . 

In particular, in this paper our ideas are tested using some very well known toy 
models of classically uniformly hyperbolic, area preserving (i.e. symplectic) maps over 
T 2 , such as the Sawtooth maps and Perturbed cat maps. Again, here we restrict ourself 
to recall the main properties needed for our considerations, refereing the reader to the 
cited references for additional mathematical details. While the first model (Sawtooth 
map) deals with piecewise linear discontinuous maps that develop diffractions effects 
during the evolution of given initial quantum states, the Perturbed cat maps are a 
prototype of non linear smooth chaotic systems and we are interested to explore the 
effects on QCF of both of these phenomena. 

For any given K G R, the classical Sawtooth map is defined on the torus 
T = [0, 1] x [0, L], L G N, with the map 4> s : (q,p) i-> (q',p'): 

p> = p + Kq (mod L) , 
q' = q + p' (mod 1) . 
|| Defined in general for any given density matrix. 
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This map (|27|) is not smooth for K ^ Z, it is ergodic and exponentially mixing 
with respect to the Lebesgue measure and it has a simple expression for the maximal 
Lyapunov exponent: 



The corresponding quantum evolution operator Us G C over the Hilbert space of 
size N is defined as [23] 



where we introduce operators m\p m ) = m\p m ) and ri|g n ) = n\q n ) with n,m e Z N . 
Please note that for odd N periodicity of free-propagation in direction of momenta 
is broken unless the dimension of the phase space in that direction, L, is even. The 
classical discontinuities of the dynamics strongly affect the propagation of states and 
the semiclassical properties of the system. Because of this, only few partial rigorous 
Egorov estimates can be given in a time regime that, as we will discuss, it is really too 
short for fully exploring our observed numerical results [23 E3] • 

We now turn to the so called Perturbed cat map over T 2 = [0, l] 2 . We start from 
the very well known generalized Cat map <fi : T 2 — > T 2 represented, in the natural 
coordinates (q,p), by the matrix 



with |Tr(0)| = \a + d\ > 2 to ensure hyperbolicity. As a consequence, <fi has two 
eigenvalues of the form e ±A , where A > is the uniform Lyapunov exponent of the map. 
On the plane, the two corresponding unstable and stable eigenspaces are generated 
by two vectors v u and v s of the form v u = (l,u) and v s = (l,s), with quadratic 
irrational slopes. If 7r : M 2 — > T 2 denotes the usual projection mod 1 on the torus, 
the image tt (R ■ v u ) coincide with the unstable manifold W u (0, 0) passing trough the 
origin. Similarly this holds also for the stable manifold W s (0,0) = it (R ■ v s ) and 
this linear structure will be immediately lost when one introduce some nonlinearity 
by perturbation. Because of the irrationality of the slope, each invariant manifold is 
dense on the torus and this property remains under perturbation. This implies that 
longer segments of W u (0, 0) becomes arbitrarily close to other pieces of itself and to 
the origin^. This phenomena clearly will produce interference effects in the long time 
quantum evolution of initially localized states and it will be responsible for the behavior 
of the observed QCF at later times. Any such linear hyperbolic map represents the 
most elementary example of two dimensional discrete highly chaotic dynamical system: 
in particular it is ergodic with respect to the Lebesgue measure, it has a dense set of 
unstable periodic orbits and it is exponentially mixing. More precisely, for any given 

% In the linear case, and then also in the perturbed one, elementary diophantine estimates show that 
two different branches of a piece of length I of W u are at distance > y, for a certain constant C > 0. 



A s = log -(K + 2 + y/K(K + 4)) 




(28) 




(29) 
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observable f,g £ C k (T 2 ) with zero average and given finite regularity fceN: 



r 2 



< Ce- tfcA || /|| c * |M| C *. 



This exponential decay of (classical) correlations and its behavior under nonlinear 
perturbations will clearly play a fundamental rule in understanding and explaining the 
numerical observed quantum-classical fidelity. 

We will now consider perturbations of this linear map by composing it with a time 
p > flow of a global Hamiltonian H(q,p) on T 2 . Assume if to be a real function of 
the torus of the form: 

H(q,p) = 22 a n m sm(27r(nq + mp)) + b nm cos(2n(nq + mp)) , 

with rapidly decaying coefficients a nm , b nm . Then H induces a vector field X H = 
— if') on ^ 2 anc ^ a corres P on ding Hamiltonian flow 6 t : T 2 — > T 2 , 9 (q,p) = (q,p). 
We now consider the time t = p > as a perturbative parameter and define a Perturbed 
cat map either as M := o # M or M := 6*^ o 0. It is known |2Sl EE] that there exists a 
/Vax > 0, which depends on and if, such that for values of fi G [0, /x max ), the map M 
is still an Anosov map of the torus and moreover it is C°-conjugate with the linear map 
0. Namely, for each p, fixed there exists a continuous Holder function ^ on T 2 such 
that ty' 1 o o vp^ = 0^. In particular this implies that both the rigid spatial structure 
and the stability exponents of the periodic orbits can change drastically, but the global 
topological entropy remains constant. 

In this paper we consider the particular perturbations for which the Hamiltonian 
function H of the perturbation depends only on one coordinate. In this case, the 
corresponding Hamiltonian flow gives rise to non-linear shears. More precisely, assume 
for example 

H (q) = a n sm(2nnq) + b n cos(27mg) . 
It is then easy to see that the corresponding Hamiltonian flow at time p is of the form: 




p ) ' \ p + vf(q) 

a shear in momentum, where f(q) = — ^p- Similarly, a Hamiltonian of the form 

H(p) = a m sin(27rmp) + b m cos(2nmp) 

meZ 

generates a shear in position 

9 1 ._ I V + Mip) 
P 



(30) 




(31) 
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with g(p) = — — . For our numerical explorations in order to simplify the numerical 
dp 

procedures without loosing generality, we consider the torus T = [0,1] x [0, L] (L EN) 
and we perform our calculations with the Perturbed cat map (fi^ : (q,p) (q',p') given 
by the formula 

p' = q + p - fi s in(27rg) (mod L) , 
q 1 = q + p' (mod 1) . 

The hyperbolicity of this map is not uniform and the Lyapunov exponent Ap(/i) changes 
with perturbation strength fi as shown in the figure (£[]) . Also in this case, the quantum 
counterpart of the classical map (}3*2~|) and some of its semiclassical and spectral properties 
have been extensively studied [23 I2E1 123 EOj- In particular the evolution operator 
Up E C NxN reads 

Up = exp (-i^m 2 ) exp (i^n 2 + i|| cos 0£n) ) , (33) 

To successfully develop the problem of correspondence and to study the observed 
behavior of the QCF we need to explore the theory of coherent states on T 2 and their 
evolution under nonlinear or discontinuous hyperbolic dynamics both in configuration 
and phase-space representation (Wigner function). As we will briefly report here and 
in more mathematical details in [31], recent important estimates have been obtained in 
[3*U] for what concern the evolution of localized coherent states under the perturbed cat 
dynamics. It is worth mentioning that even if these represent the best known estimates 
about time propagation of states, they have a genuine semiclassical feature and they 
rely on a detailed control on the difference between the classical and quantum evolution 
(Egorov Theorem). For these reasons they do apply at a time scale where both non 
linear classical effects and quantum interference are negligible, allowing us to rigorously 
prove only the observed non decaying of the QCF for short initial time. This of course 
excludes the possibility of understanding the "interesting" part of the QCF behavior 
at larger time, where different mathematical approaches are needed, as we will here 
heuristically discuss. More precisely, when translated at the level of the Wigner function, 
these rigorous "Egorov estimates" allow to gain the explanation of the behavior of the 
QCF only up to time Ti = Tg/2, or Ti = Te/6, depending on whether the map is linear 
(with possible discontinuities), or non-linear, respectively. Moreover, the breakdown of 
these estimates also offer some understanding of the mechanism which produces the 
appearance (at least heuristically) of the observed exponential decay controlled by the 
maximum Lyapunov exponent, up to the 1/N plateau at time T2 = Ti + Te, as we will 
discuss more deeply at the end of the section. 

Let us here briefly recall the main definitions and estimates concerning coherent 
states and their evolution, refereing to [30J for the technical mathematical details. 
As a starting point we consider a localized coherent states at the origin denoted by 
^o,r G L 2 {R), defined as [32] 



*o,r(g) 



\\ I/ 4 



Txh~ 
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Figure 1. Lyapunov exponent of the classical Perturbed cat map as a function of the 
perturbation strength, obtained as an average over 100 trajectories of length 10 6 with 
random independent initial points. 

Here r G C is a point in the right half-plane related to the squeezing of the coherent 
state. Coherent states ^ X:(T localized on other points x = (q,p) in phase space can then 
be obtained by the action of the Heisenberg translation operators. The corresponding 
coherent state ^f X)T j2 G TLn on T 2 is then obtained by projecting down the state through 
the operator P T 2 as defined in the eq. (J22j) |2~U] : 

To be more concrete, if r = io~ with width o > then the coherent state at the point 
x = (qo,Po) G R 2 , denoted again with a small abuse of notations by ^ Xja j 2 , becomes for 
h= 1/(2ttN): 

^ (£) = EM^ +z/ ) 

1/4 

e i27mp ^ e -^(n-q N+uN) 2 e i2irNp u 

Note that if go ^ 1 — then in the classical limit we can approximate the torus 
coherent state by the planar one 20J: 

In the study of the quantum-classical correspondence we need to compare the quantum 
and classical evolution of the classical probability density corresponding to the initial 
phase-space representation of a given coherent state. As we will discuss also later 
on, even if there is some freedom in the choice of the Wigner representation on the 
torus, the previous estimate yields that the classical probability density that closely 
resembles the phase space representation of the coherent state ^ X:(7 j2 turns out to be 
a periodicised Gaussian function on T 2 . This last statement follows, roughly speaking, 
from the following two facts |2()j : 
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(i) The planar coherent states, given in terms of Gaussian wavefunctions are covariant 
with respect the Fourier transform F^, namely Ff^ xa = ^fx,i/oi where Fx = 
ipo, —qo)- Now when {2t[K)^ 1 = N, then P^Fn = F^Pr 2 , where Fn again denotes 
the discrete Fourier transform. This shows that the above covariance still holds also 
after projection on the torus: Fn^x,o-j 2 = ^Fv 1 ,? 2 - 

(ii) The phase space representation of the coherent state, namely its (discrete) Wigner 
function, is essentially the Fourier transform of the overlaps between the coherent 
states and its translates on T 2 , which turns out to be related to overlap between 
planar coherent states, i.e. Gaussian integrals: 

mez 2 

_ (_^JVmi"i2+mAk e 27ri(a;A(k+mAr)) e -^Q (T (m+k/JV) 

mez 2 

where m A k = mi^ — m2^i and Q CT (k) = ok\ + o~ x k\. As a consequences of this, it 
can be seen that up to exponentially small errors in N (at least around the point where 
the coherent state is localized), we can approximate the genuine initial Wigner function 
of our coherent states essentially with the periodicised Gaussian function on T 2 

/W)(?,P) = Dn f^e-^^A ^Te- 2 ™- 1 ^-^ 2 

The scalar factor D N is pinned down by normalization J T2 dg dp P( qo , Po )(q,p) = 1, which 
is in the case of continuous torus T 2 equal to Djv = 2N. 

Let us now make a non technical analysis of the semiclassical behavior of the QCF 
for the short time regime t < T±, where no decay is observed. As shown in [HO], f° r short 
time (say t < Tg/6) the evolution of ^/^^.t 2 is governed by the linearized dynamics 
around the classical trajectory of x. After that time, the Gaussian wavepacket starts 
to be seriously effected by the nonlinearity of the dynamics. To be more precise, let us 
for simplicity consider an initial coherent state ^f X)T with arbitrary squeezing parameter 
r E C and localized around the origin (fixed point of the (un)perturbed dynamics). Let 

*-(iS:)' 

be the linear hyperbolic map of the plane obtained by considering the tangent map 
of M at the origin. If we consider this as a global map on R 2 , we denote by 
the corresponding quantum unitary propagator (metaplectic representation) acting on 
L 2 (R) is written as [HI 

wlvM = 7sbsr xp te (<w _ 2q,q + a " q2) ) ' 

Now, for "short time" t, the evolution of the coherent state through the true perturbed 
dynamics Ut can then be approximated by an Ansatz := V^o tT = e iet ^o tTt , where 
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the phase 9 t and the "squeezing" r t of the new coherent states are given by the Maslov 
multiplier and the homographic transformation respectively: 

t = -~axg(a('> + rbf) , r t = , where (<%)< = ^ . 

z ay + rby \ Cn j 

Using the "projector" P T i both states and the propagator can be projected to the torus: 
V^o,t,t 2 '■= Pf 2 Vfto,T, namely ^tj 2 = -Pt 2 ^*, and the following can be proven |3Uj : 

|| E^T,^ ~ % T 2 \\ Hn < CfiW^, 

where A M denotes the maximal Lyapunov exponent for the nonlinear Perturbed cat 
map. This estimate follows from the localization condition fiAq(t) 3 <C h, where 
Aq(t) ~ h 1 / 2 e tx ^ is the width of the evolved coherent state. Below this break time 
T « the state is still very localized at the origin, being supported in a diameter 

of order ~ h 1 ^ 3 , with no interference appearing in the phase-space representation. In 
the case of piecewise linear hyperbolic maps instead, where nonlinearity is absent and 
only diffraction effects at the singularities take place, techniques similar to the ones 
developed for the Baker map [2U] should allows a precise Egorov estimates for longer 
times, up to t ~ Te/2 + . At this point serious interference effects between different 
branches of the elongate states takes place, semiclassical approximation break down and 
genuine quantum phenomena start to take place, as the developing of negative values 
in the Wigner function that are responsible for the observed time scales T 1; T 2 and the 
corresponding decay of the QCF as discussed in the next section. Without entering the 
details here, it is important to stress that, when these estimates are translated into a 
phase-space representation through the Wigner function, then one can easily infer that 
G(t) ~ 1 for times below half of the Ehrenfest time for piecewise linear maps and for 
times below pa Tg/6 for the nonlinear smooth dynamics. This result is not surprising, 
at least from the physical point of view, and we would like now to develop some better 
understanding (at least heuristically) of the "long time" behavior of the QCF F(t). In 
order to do this, we turn now to discuss the Wigner function(s) we have been using in 
our numerical investigations. 

We must stress at this point that while a rigorous and unique mathematical 
definition of the Weyl- Wigner formalism is now available (see e.g. [22]), the finite 
dimensionality of the space of quantum states and the remark following eq. (|25jl yield 
quite counter intuitive (at least from a physical point of view) Wigner function defined 
only on a 2N x 2N discrete lattice of T 2 . Even if this Wigner function does satisfies 
all the natural properties of a correct phase space distribution, such as for example eq. 
(J2HJ), and even if it gives the expected marginals or expectation values when integrated 
on linear subspaces of the phase space, it is sometimes better to use an alternative 

+ For the nonlinear case at times t > Te/6 bigger than the breaking time of the linear evolution, 
the states starts to be quite stretched along the unstable direction, and only suitable WKB Ansatz 
approximation allows one a certain control of its time evolution up to about t w Ie/2 [23] 
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version of the Wigner function which offer a more intuitive interpretation of its phase 
space overlap with the evolved Gaussian classical distribution pt(q,p)- Even if the final 
numerical results do not depend on this choice, for the sake of clarity in the presentation 
and also for the sake of generality in the interpretation of the mechanisms yielding the 
exponential decay of the QCF, we discuss here the two different approaches: the discrete 
one ( [131221) and the continuous one [T7] . As we will see, we can benefit from studying 
both formulations of the Wigner function. 

3.2. The continuous Weyl- Wigner formalism 

The Weyl- Wigner formalism introduced by Agam&Brenner ^7] defines the Wigner 
function by the following formula 

{N-l)/2 

W4n,m) = - e- i % n ' m ~6(2l-2n + n')(q l+n ,\ij)(ij\q l ), (35) 



N 

n',l=-(N-l)/2 



(n, m) £ T? N , 

with the physical position on the torus T 2 = [0, l] 2 given as (Jr, ^) e T 2 and where 

S(k) = - V = 1 sin(7rfc/2) 

N ^ Nsin(Trk/2N) ' 

m'=-(N-l)/2 K 1 ' 

This formulation of the Wigner function requests that the Hilbert space dimension N 
is odd. This could be a deficit, as we have no concrete reason for it, except maybe that 
it is rather natural to assume the same number of negative as positive wave-lengths 
resulting in the odd number of all wave-lengths. The Wigner function given by (|36|) has 
the following properties 

N-l N-l N-l 

J2W4n,m) = \(p m \ij)\ 2 , J2w^(n,m) = \(q n \ij)\ 2 , ^W^(n,m)=l. 

n=0 m=0 n,m=0 

The full Weyl- Wigner formalism that is incorporated in jTTj can be presented by using 
the kernel/point operator A nm for projecting operators on the function over phase-space 
grid and its reconstruction from it 

N-l 

A = -A ^ a nm A nm , a nm = tr{j4 nm y4.} , 

n,m=0 

(N-l)/2 

A nm = - Yl e " ¥n ' m ~ 6 ^ 1 - n ) + n ') l»> • 

n',l=-(N-l)/2 

This type of Wigner function can be extended to the continuous torus (q, p) G T 2 by 
the following definition 

W4q } p) = W^(Np,Nq) , 
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with the basic properties 

J ^ dqdpW^(q,p) = , J ^ dq dp W${q, p) = , 

which can be treated as normalization of the Wigner function. We only use presented 
Wigner function in the continuous formulation as this is a more compatible object with 
the classical density. The classical-quantum fidelity of the evolving Wigner function W^t 
and classical density p* = p°((p~ t (q,p)) on the continuous torus is defined as 

F(t)= [ dqdpW^(q,p)p t (q,p), (36) 

= [ dqdpW^(4>\q,p))p%q,p), (37) 

where we use the phase space conservation of the map 0*. Here we obtain an interesting 
object to discuss H^/,t(0'(g, p)) that we call the classically echoed Wigner function and 
corresponds to z l (x) discussed in the Euclidean geometry. By using this object we can 
directly observe decay of classical-quantum correspondence on the phase space as it 
spreads from original Gaussian form over the whole phase space. This is nicely depicted 
in figure El where we plot W^t(q,p), p l (q,p) and W^t^^q^p)) parallel to each other 
at each time step in the saw-tooth systems. On figures representing classically echoed 
Wigner function W^t^^q^p)) we can clearly see how is this function stretched in stable 
directions of the classical map as predicted by the analysis on Euclidean phase space. 

We numerically calculate the quantum-classical fidelity F(t) or its relative value 
G{t) = F(t)/F(0) for the Sawtooth map system and Perturbed cat map. The results 
for different parameters are presented in the figures El for the Sawtooth map and in the 
figure E]for perturbed cat map. The results are averaged over different positions of initial 
Gaussian distributions to get rid of features depending on individual initial states. In all 
results we can clearly see the exponential decay of relative QCF (G(t)) ~ e ~ Amax< with 
the rate equal to the maximal Lyapunov exponent A max down to the ergodic plateau 
(G(t)) c = jr- An exception is the case of the Sawtooth system at small parameters K, 
where the expected break time t^ is rather big and we can also experience the correlation 
decay. Obtained result are all in the frame of presented theory for the Euclidean phase 
space. Note that in order to obtain reliable numerical results we have to use more then 
TV 2 points in integration over the phase space, where N is the Hilbert space dimension. 

3.3. Discrete Weyl- Wigner formalism 

The other Weyl- Wigner formalism that we use is described by Miquel et.al.jTHj, by which 
the Wigner function reads 

W4n,m) = ^—Y^^\qn-k)(q k \^)e-^ km , (n,m) G Z 2 2N , (38) 

k=0 

where the physical position on the torus T = [0, l] 2 is given with (^, ^) e T 2 . Even if 
being defined on a 2N x 2N lattice, turns out to be completely determined by iV 2 
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-0.001 -0.0005 0.0005 0.001 0.0015 0.002 



Figure 2. The density plot of the Wigner function (first column), the classical density 
(second column) and the classically echoed Wigner function (last column) at t-th step in 
the Sawtooth map systems at parameters K — 0.5, L = 2 and Hubert space dimension 
N = 101. The white lines indicate unstable and stable direction of the symplectic flow. 

of its values, as the following symmetry shows: 

W^n + a q N,m + a p N) = W^{n,m){-iy" m+ ^ n+a ^ N 
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Figure 3. The average QCF for the continuous Wigner function formulation in the 
sawtooth map at N = 101, 201, 401, L — 2 and different K = 0.1, 0.2, 0.5, 1. The 
average is calculated over 10 different initial Gaussian coherent states in the case 
N = 201, 401 and over 50 in the case N = 101, respectively. The fitted parameters are 
ui = 0.978 and v 2 = 0.935. 




Figure 4. The average QCF for the continuous Wigner function formulation in the 
Perturbed cat map for AT=101, 201, 401, L = 2 and different fj, = 0.1, 0.5 (left,right). 
The Lyapunov exponents here in use are Ap(0.1) = 0.964 and Ap(0.5) = 0.952. The 
averaging is done as for results in the figure 



with ci qp G {0, 1}. Moreover, this Wigner function satisfies all the main natural 
requirements for a correct phase-space representation. In particular it does have the 
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right "Lagrangian averages" on the points in phase space on the lattice of size N: 

2N-1 2N-1 2N-1 

Y^W4n,m) = 5 2 J:-\(p k ^)\ 2 , W,p(n,m)=5 2 n k -\(q k n\^)\ 2 , ^ W^(n,m) = l 

n=0 m=0 n,m=0 

The Weyl-Wigner formalism that produces this Wigner function is generated by a point 
operator with a simple algebraic form 

2JV-1 

^n,m tr{-A-A n)TO } , A N / " ii .in ■' ' ii .! ii ■ 



ti{AA n ^ m } , A — N ^ ^ a n , m A ri 

n,m=0 



fc=0 

In the presented formulation we perform numerical experiments considering the discrete 
version of the quantum-classical fidelity, 

2JV-1 
n,m=0 

where we sum the overlap of the Wigner function and classical density over the lattice 
2N x 2N. This formula for fidelity ()H9j) has the time complexity of 0(N S log(iV)) and is 
by a factor O(N) smaller than the one using the Wigner function defined by Agam & 
Brenner. This enables testing of classical-quantum correspondence and its break-down 
at higher dimensions. We have studied such QCF in different systems: perturbed cat- 
map, saw-tooth map and also baker map and general linear skew shift (not reported 
here) 

N-l N-1 N-1 

J2w4n,m) = \( Pm \^)\ 2 , ^^(fi,m) = |(? B W| 2 , ^)W^(n,m)=l. 

n=0 m=0 n,m=0 

The results are shown in the figure El for the Sawtooth map and in the figure |H1 for 
Perturbed cat map. Generally, in the generic case where the lattice is not conserved under 
the classical map, and the trajectories starting at lattice points are scattered over the 
whole phase space, we get the results for fidelity which look asymptotically (as N — > oo) 
the same as in case of the Agam & Brenner's definition of Wigner function. There exist 
also non-generic interesting cases. For instance in the Baker's map at dimension iV = 2 ra , 
where the classical map in time flattens the lattice onto the ordinate axis. The other 
example is the saw-tooth map at K 6 N called the cat map, which conserves the given 
lattice and the fidelity is constant, meaning it does not decay. This is a consequence of 
the fact that the cat map is Egorov exact [TT?] . 

The discrete definition of the Wigner function enables to study in more detail 
characteristic times of the QCF time dependence. We can clearly recognize two distinct 
time scales. The first one, denoted by T 1; is the time at which F(t) starts to strongly 
deviate from initial value F(0), entering the regime of exponential decay. The second 
time scale, denoted by T 2 , is the time when F(t) exits the exponential decay and starts 
to relax towards the plateau. Actually, we define the times T% and T 2 using the average 
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Figure 5. Time evolution of the average QCF for the discrete Wigner function 
formulation on the lattice 2N x 2N in the Sawtooth map at L = 1 for K=0.1, 0.2, 0.9, 
1.1. The fitted paramater /ii = 0.945. The average is taken over 500 initial Gaussian 
coherent states. 




Figure 6. The average QCF for the discrete Wigner function formulation in the 
Perturbed cat map at L = 1 and different values of the parameter fi = 0.1, 0.5 (left, 
right). Lyapunov exponents are the as used in figure 0] The average is calculated over 
100 different initial Gaussian coherent states. 



relative QCF, averaged over different initial states (G(t)) and study the exponential 
slope of (G(t)): a very good agreement with the exponential fit 67 fit (t) = exp(7\ — A max t) 
is then observed for large N, where A max is the maximal Lyapunov exponent of the 
classical map. We consequently define the times T\ and T 2 as the times when the fitting 
function G^ t {t) hits the initial value G(0) = 1 and the plateau 1/N respectively as shown 
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Figure 7. Scheme for precise definition of two time scales, T\ and T2, of decay of 
relative QCF G{t). The thick straight line is the best exponential fit in the region of 
decay. 



We see that the times T 2 and Ti are connected by the relation 

X max (T2-T 1 ) = \ogN (40) 

Because of that trivial relation we discuss only the time scale T\ in the following. In the 
simulations done with the Sawtooth map and Perturbed cat map, see figure El we found 
that T\ is approximately described by 

A max T\ « AlogN + B (41) 

where A and B are constants, depending on the model but not on N. We are are 
interested only in leading term and the coefficient A. By numerical analysis of the 
models at different values of parameters we found that A in Sawtooth map T2 is almost 
independent of A max and is near to A = 0.5 ± 0.05, but in the Perturbed cat map 
the coefficient A slightly varies with the coefficient \x and is approximately equal to 
A = 0.2 ±0.05. 

Another interesting aspect of dynamics on the torus is that we can clearly see the 
connection between decay of QCF and the negativity of Wigner functions. In the figure 
E]we show some numerical results where we plot together the relative QCF G(t) and the 
relative fraction of phase space P_ (t) supporting the Wigner function with the negative 
value, starting from the same initial Gaussian packets. From a detailed numerical study 
of statistical properties of the Wigner function for chaotic maps on compact phase spaces 
[T4*] we know that P-(t) should saturate in time to the value predicted by the random 
wave model 

P- rand = - - -L T C^^dx = - - -£= + 0(r 2 ) (42) 
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Figure 8. Dependence of the time Ti in decay of average relative QCF (G(t)) on 
the Hilbert space dimension TV for the Sawooth map (a) and Perturbed cat map at 
parameter values which are indicated in the figure. Inserted solid lines are linear fits: 
-0.234 + 0.508 log TV, 0.995 + 0.500 log TV, 1.411 + 0.519 log TV in the case (a) and 
2.539 + 0.159 log TV, 1.944 + 1.169 log TV in the case (b). The average in (G(t)) is taken 
over 500 initial Gaussian packets. 



where r = W/a w is the ratio between the average Wigner function W and its standard 
deviation a w . For the considered Wigner function formulations on the torus we have 

continuous : r = (N - 1)~ 1/2 , discrete : r = (4tV - 1)~ 1/2 . (43) 

We see that decrease of QCF closely follows the increase of P-(t), but there is obviously 
a strong difference between the discontinuous piece-wise linear (i.e. Sawtooth map) 
and non-linear smooth maps (i.e. Perturbed cat map). Roughly speaking, the QCF 
decay appears when the Wigner function starts to develop, due to genuine quantum 
interference effects, Gaussian statistics for its values distribution. This correspondence is 
sharp for linear maps with discontinuities, like the saw-tooth map, whereas for nonlinear 
continuous maps, like perturbed cat map, the relaxation of Wigner function statistics 
happens little later than the time scale Ti of QCF decay. Note that such behaviour has 
been observed for several other, different values of TV and as well as different values of 
parameters of the two representative models, not shown in the figure. 

4. Conclusion 

In the present paper we have discussed global phase space aspects of quantum-classical 
correspondence. We have defined a quantity, called quantum-classical fidelity (QCF), 
which characterizes L 2 norm distance between time evolving classical and quantum 
phase space distributions, starting from the same initial function, e.g. Gaussian. 
Quantum phase space distributions are described using Wigner- Weyl formalism. 

Using heuristic arguments and numerical simulations, we conclude that QCF is 
characterized by two Ehrenfest-type time scales for classically chaotic dynamics and 
effectively compact phase spaces. QCF remains close to its maximal value up to time 
Ti, and after that starts to decay exponentially until it saturates to a plateau given by 
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Figure 9. Time evolution of the average relative QCF (G(t)) and of the average 
portion of negative WF at N = 2048 in the Sawtooth map at K = 0.9 (a) and 

Perturbed cat map at /i = 0.1 (b). The averages are taken over 500 initial Gaussian 
packets. 



finite effective Hilbert space dimension, at time T 2 . Both time scales are proportional 
to logarithm of effective Planck constant. Our results strongly suggest that the rate of 
exponential decay of QCF, after Ti, is given by the maximal Lyapunov exponent of the 
corresponding classical dynamics. 

We have given heuristic arguments for our results in the case of Euclidean phase 
space geometry, and in addition, we have provided numerical results for generic chaotic 
maps on compact toroidal phase space. We outline subtle but significant differences 
between linear discontinuous maps and non-linear maps. 

We have suggested that understanding of QCF may provide a key to the 
correspondence between quantum and classical Loschmidt echoes, namely we have 
proven an inequality which shows that quantum Loschmidt echo simply follows the 
classical one up to time scale T\. 

Several open problems and directions of future work can be suggested. Rigorous 
work is needed to understand the above results. One possible direction is the perturbative 
treatment of the Egorov property, like for example in the perturbed cat map studied 
above. Further, one may investigate the behavior of QCF for integrable systems and 
(weakly) chaotic systems with more complicated classical phase space structure, for 
example of non-uniformly hyperbolic or even KAM type. 
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